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Abstract. We consider third order Lovelock gravity coupled to an U(l) gauge field for which its Lagrangian 
is given by a power of Maxwell invariant. In this paper, we present a class of horizon flat rotating black 
branes and investigate their geometrical properties and the effect of nonlinearity on the solutions. We use 
some known formulas and methods to calculate thermodynamic and conserved quantities. Finally, we check 
the satisfaction of the first law of thermodynamics. 

PACS. PACS-key discribing text of that key - PACS-key discribing text of that key 


1 Introduction 

Lovelock gravity is one of the higher derivative gravity 
theories, natural generalization of Einstein’s General Rel¬ 
ativity, introduced by David Lovelock ^\m\ in 1971 

. d/2 

la = / (1) 

d k=0 

where (a^) is an arbitrary constant and {£k) is the Euler 
density of a 2fc-dimensional manifold 

P, — XklVl - PkVk DPlO-l UPkCTk ('9'\ 

In this equation, is the generalized totally anti¬ 

symmetric Kronecker delta and {RfZ) is th® Riemann ten¬ 
sor. 

It has been shown that gravitational field equation, 
arising from variation of the action [U only contains the 
terms with at most second order derivatives of metric 
which indicates the quantization of linearized Lovelock 
theory is free of ghosts [IS]. Also, such theory may be 
used in the context of AdS/CFT correspondence to inves¬ 
tigate the effects of including higher-curvature terms jS] 
misiig) . Another interesting property of this theory comes 
from the fact that there is a delicate relationship between 
Lovelock gravity and low energy effective action of het¬ 
erotic string theory in addition to the Calabi-Yau com- 
pactihcations of M-theory [ills] . 

On the other hand, it has been shown that loop correc¬ 
tions of gravity [lUKIIllIllIllIi] may be lead to a nonlin¬ 
ear Born-Infeld type Lagrangian [TSlfTS] . Hoffmann was 


the first one to attempt to couple the nonlinear electro¬ 
dynamics with gravity m- One of the current interesting 
subjects in black hole physics is the investigation of the 
effects of nonlinear electrodynamic fields on the gravita¬ 
tional field 

[SSlISillSSlIdbllST] . Considering nonlinear model of electro¬ 
dynamic fields provide a powerful laboratory to analyze 
the properties of the black hole solutions. 

In this paper, we obtain d-dimensional rotating black 
brane solutions of the first four terms of Lovelock theory 
which are cosmological constant, Einstein, Gauss-Bonnet 
and third order Lovelock terms, in the presence of nonlin¬ 
ear electromagnetic field and analyze their properties. We 
should note that third order Lovelock term is an Euler 
density in six dimensions and in order to have the con¬ 
tribution of all the above terms in the field equation, the 
dimension of the spacetime should be equal to or larger 
than seven. These solutions can be regarded as a gen¬ 
eralization of the Lovelock-Maxwell gravity [551155] . We 
discuss some geometrical and thermodynamic properties, 
like curvature, temperature and entropy. We also make 
some comments on the effects of nonlinearity parameter 
and Lovelock coefficients. 


2 Field Equations 

The action of third order Lovelock gravity in the presence 
of power-Maxwell-invariant source may be written as 

T = —— f d'^xy/ — g\R — 2A -|- (X2C-2 T 0:3/13 -|- 
IStt Jm 

(s:.F)"]+Ib, 
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where R is the Ricci scalar, A is the cosmological constant, 
Oi’s are Lovelock coefficients, £2 and £3 denote the Gauss- 
Bonnet Lagrangian and the third order Lovelock term, 
given as 


£2 = (4) 


_ O D DPT I Q DP^ D<TK DPT 1 

3 - ap^ p.n + 

r\A TDUUCrKs TD D/O 1 Q D TDUU<TK TD 1 O D D DP^ 

z4xt lL(ji^i/pI\ ^ “T OlLlt iZlXlXpi/li, 

. (5) 

The last term in the first integral is the Lagrangian of 
power-Maxwell-invariant theory, where k is an arbitrary 
constant and R is the Maxwell invariant which is equal 
to F^^ = d^Av — duA^ is the electromagnetic 

tensor held and is the vector potential) and the last 
term is the boundary term which is chosen such that the 
variational principle is well dehned. Varying the action © 
with respect to and one can obtain the gravita¬ 
tional as well as the electromagnetic held equations 


-I- — 


pLV 

-2 


pLlA 


KsF^pFS {kFY-^ - {kFY 


( 6 ) 


^^.(y^g{nFr-^F^^-)=Q, (7) 

where = Riiv ~ ^9^^ is just the Einstein tensor, and 

(2) (3) 

G/tj; and GfiJ are the second and third order Lovelock 
tensors obtained as 


G^fj = 2{R^,^rRr^ 




‘^‘Rfj.pvaR'’'^ ~ 2,RpaR'^i, + RRfw) 

( 8 ) 


G piy of4_/l R(TKXpR i/rp 8ri X(jR rpR UpK ' 

9 P TaK TD D^P _ TDTpaK. TD TD 1 

ZJt^ rLfjf^XpiL IT ^aKTp^ypi “r 

O DT D^f^ DP I Q D^ DTp DK 
l/(7p^ T/i-^ K ' UTK^ api^ p 

+mj^^R,^ppRY - mj^^R„^rpR% + 
m^P’^^R.^rpR.p + 2RR^^^PRrp.p + 
8R\^^RP,R% + 8R\^pRP,R% - 
8R%,pR\RP^ - 8RZpR%R.p - iRRYppRY 
-\-4:R^pRptR,^p — 8R^jjRrpRPp A 4 lRR^pRP p 

—RYRvp) — -C^gpu- (9) 


£36 = 3(P - 2Gi"^iL“'’ - 12P,bJ“'' + 2RJ - 

AKRabcdK^^K^^' - 8 RabcdK‘^^KlK^Y- (12) 

In these equations jpi, and K are, respectively, induced 
metric and the trace of extrinsic curvature of boundary, 
G^^J and G^^jj denote the {d— l)-dimensional Einstein and 
second order Lovelock tensors (Eq. ([5|)) of the metric jab 
while J and P are the traces of 

Jab = ^{2KKacKt + KadK^‘^Kab - 

2KacK^‘^Kdb-K^Kab), (13) 

and 

Pab = 1{[K^ - 6K^K^‘^Kad + 8KKadKiK<^^ - 
QKadK'^-KafKf- + 8{KadK-Y‘"]Kab - 
[AK^ - 12KKedK^‘^ + 8KdeK}KfYKacK^ - 
2AKKacK^'^KdeKl; + 12{K^ - KafK^YKacK^^Kdb 
+2AKacK^‘^KdeK^fKtf}. (14) 

3 The d-dimensional Charged Rotating Black 
Branes 

Now, we should consider a rotating spacetime and investi¬ 
gate its properties. Since the rotation group and the num¬ 
ber of independent rotation parameters in d-dimensions 
are, respectively, SO{d — 1) and the integer part of (d — 
l)/2, one of the rotating metrics with fc < [(d — l)/2] ro¬ 
tation parameters and zero curvature boundary may be 
written as 

/ V dr^ 

ds"^ =-f{r) i^dt-X^aid(j}ij + + 

9 /c 9/2 

^ ^ ((3-2(it i^l d(l)i^ ^ cZjfi^^) 

i—l i<j 

d—2 — k 

+r^ Y, dxl (15) 

2=1 


where n = + ^'^‘1 Hi® angular coordinates 

are in the range 0 < (pi < 2 tt. Using the suitable gauge 
potential ansatz 

Ap = h{r) {F5^ — aiSY) (no sum on i), (16) 


The Lovelock supplementation, Ib, is Gibbons-Hawking- ^ following differential equations 

York boundary term which may be written as 

^ . Eg?, = (2s - l)rh"-b (d - 2)/i'= 0 , (17) 

Tb = --^ -b a2£2h + Q!3£3b], (10) 

JdM with the following solutions 


where 

£2b = 2(j-2GWE“''), 


( 11 ) 


h{r) 


—glnr for s = 

qp(2s-d+i)/{2s-i) otherwise ’ 


(18) 
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where prime and double prime denote first and second 
derivatives with respect to r, respectively and q is an in¬ 
tegration constant. It is worthwhile to note that, one can 
choose s = d/A to obtain inverse-square electric field. In 
this case the expression of the electric field does not de¬ 
pend on the dimension and its value coincides with the 
Reissner-Nordstrdm solution in four dimensions. In ad¬ 
dition, in arbitrary dimensions with s = {d — I)/2, the 
gauge potential is logarithmic and it coincides with the 
charged BTZ solution. In order to have a sensible asymp¬ 
totic structure, the electromagnetic field should vanish for 
large values of r. This condition leads to s > 1/2. In this 
paper we would like to investigate the solutions of general 
case (i.e. s > 1/2 and s ^ {d— l)/2). In addition, one can 

find that T = —2 ^ and the power of Maxwell in¬ 

variant, {kjFY , may be imaginary for positive k, when s is 
fractional. Therefore, we set k = — 1 to have real solutions 
without loss of generality. 

To find the metric function /(r), we should consider 
the components of Eq. (jH]). At first, we consider the static 
case {at = 0). One may show that the tt-component of Eq. 

may be written as 


Eqtt 


d-2 




+ [6!(V)«3/"-4!a2C 


:-2 

4 


+Ar^ + 


(2s- 1) (2/i'2)V6 


yf + 2(d-2)r^] f 

yf + 2(‘^-yyl 

■ = 0, (19) 


where (y = Other nonzero components of Eq. 

® are a combination of , Eqtt and Eqh ■ For example 
a;ia:i-component of Eq. ® can be written as 


E(JxiXi — 7 

dr 


d ^ yd-8) ^ sr^ (2h'y 


-Eqtt - 


h' 


Eqh = 0. 


Since Eqh vanishes for obtained /i(r), it is sufficient to 
solve Eqtt = 0. For simplicity, we choose a special case 


as = 


(d — 3)(d — A)al 


3(d — 5)(d — 6 ) 

The solution of Eq. m can be written as 


( 20 ) 


fy = 


(d — 3)(d — 4)^2 - 


l-g(r)i/3 


where 


( 21 ) 


^ ^ ^ ^ 6(d — 3){d — 4)^2^ ^ 3(d — 3)(d — 4 )q;2M 


L{r) = 


{d-l)(d-2) 

3(c? — 3)(d — 4 )q!2 (2s — 1)^ L(r) 
(d-2) {2s-d+l) ^ 

2q^ (2s — d + 1)^ \ 

(2s - l)2r2(d-2)/(2s-l) I > 




( 22 ) 


and M is an integration constant. We should note that, in 
addition to static components of Eq. dH]), Eq. (1^ satis¬ 
fies all nonstatic components of Eq. (jH]). In addition, it is 


notable that this solution reduces to the solution of Ref. 
[39] for s = 1, as it should be. 

Due to the fact that the power of r in the denominator 
of L(r) is bigger than unity for arbitrary d > 6 and s > 
1/2, L(r) vanishes for large values of r. Therefore in order 
to investigate the asymptotic behavior of the solutions, 
one can consider the vacuum solutions (M = q = 0) 


ly 



-i I 6(cZ—3)(ti—4)a2^ 
I + id-l)(d-2) 


(d - 3)(d - 4)a2 


(23) 


Considering the former equation, it is easy to find that the 
asymptotic behavior of the solutions is AdS for negative 
A. One can put A = —(d — l)(d — 2)/2P in Eq. to 
find the effective cosmological constant of the solutions 


(d-l)(d-2) 1 - 


Aeff = — - 


1 - 


3((i—3)(d—4 )q:2 
P 


2(d- 3)(d- 4)a2 


(24) 

Now, we look for the essential singularity(ies). It is easy 
to show that the Kretschmann scalar of the metric ca) is 


= /"' + 2(d-2)^+2(d-2)(d-3)^. (25) 

After some algebraic manipulation, one can show that 
the Kretschmann scalar (na with metric function (EU 
diverges at r = 0 and therefore there is an essential time¬ 
like singularity located at r = 0. In order to investigate 
the existence of horizon(s), one should find the root(s) of 
y = /(^) = 0. In addition, we can use the fact that the 
temperature of the extreme black brane is zero. In the 
case of extreme black brane solution, the only real root of 
f(r) is its minimum. So, it is easy to find that the mass 
parameter for this case becomes 


(d_l)(2s-l) 

l2 = (d-2)] 


4As 


-2A 


^ext — 


(2i,-l)2 


(2s-l) 


(2s — d + l)(d — 1) 


(26) 


Finally, we should note that obtained solutions may be 
interpreted as black branes with inner and outer event 
horizons provided m > mext, extreme black brane for 
m = TOext and naked singularity otherwise (see Fig. [Tjfor 
more details). In addition, it is desirable to investigate the 
effects of the Lovelock parameter, a 2 and the nonlinear¬ 
ity parameter, s. For example, figures [5] and [3] show that 
both a 2 and s affect on the minimum value of f(r). Also, 
these figures confirm that in spite of a 2 , the nonlinearity 
parameter, s has a considerable role on the values of inner 
and outer horizons. 


4 Thermodynamic and Conserved Quantities 
and the first law 

Now, we are in a position to calculate the thermodynamic 
and conserved quantities of the solutions. The Hawking 




























4 


S. H. Hendi et al.: Third order Lovelock black branes in the presence of a nonlinear electromagnetic field 



Fig. 1. f{r) versus r for d = 10, 5 = 1, d = —1, 02 = 0.1, 
s = 2, and M = 3.52 < Mext (bold line), M = 4.52 = Mext 
(solid line) and M = 5.52 > Mext (dashed line). 



Fig. 3. /(r) versus r for d — 10, q = 1, A = —1, M = 5, 
02 = 0.1, and s = 2 (bold line), s = 2.1 (solid line), s = 2.2 
(dashed line) and s = 2.3 (dotted line). 



(28) 

Using the fact that the temperature of the extreme black 
branes vanishes, one can write the horizon radius of the 
extreme black brane 

2s{d-2)/{2s-l) ^ - g) ( 2 gLt (2s -d+lf \ 

A [ (2s-1)2 J ■ ^ ^ 

In order to obtain the electric charge of the solution, 
we use the generalized Maxwell equation ©• Calculating 
the flux of the electric field at infinity leads to the electric 
charge per unit volume Vd -2 


Fig. 2. /(r) versus r for d = 7, q = 1, A = —0.1, M = 1, 
s = 2, and 02 = 0.05 (bold line), 02 = 0.2 (solid line), 02 = 0.5 
(dashed line) and a 2 — 2 (dotted line). 


The electric potential U, measured at infinity with respect 
to the horizon, is calculated by the following relation 

U=A^Xnr^o.-A^X'^l^r+, (31) 


temperature and angular velocities of the event horizon of 
the black branes may be calculated by requiring the ab¬ 
sence of conical singularity at the horizon in the Euclidean 
sector of the solutions. Considering the analytic continu¬ 
ation of the metric by setting t —> it and at —> iai, 
and regularity at event horizon (r = r_|_) help us to obtain 
the Hawking temperature and the angular velocities of the 
black branes 


T+ = 


f{r+) 

47r,^ 


r+ 


(1 - 2s) 


2q^(2s-d+l)^ \ 

4(d — 2)7rS 


- 2A 


where X = + is the Killing vector. The electric 

potential will be 


jj ^ ^^Cs-d-H)/(2s-l) 


(32) 


Now, we should calculate the total mass density of the 
black branes. We may obtain finite mass through the use 
of the counterterm method [40]. One can show that for 
obtained solutions with flat boundary, Rabcd{l) = 0, the 
finite action is 


^finite — 'I' 


jd-l 


Stt , 


“7 


dM 


d-2 


,(27) 


(33) 
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It is notable that that this counterterm has exactly the 
same form as the counterterm in Einstein gravity for a 
spacetime with zero curvature boundary in which I is re¬ 
placed by leff- After some algebraic manipulation, one 
can show that 4// is given by 

^ ^ 15v^(l-A)(d-3)(d-4)a2 

3(d-3)(d-4)a2^^^^ 

P ) ’ 

which reduces to Z as a goes to zero. Using Eqs. and 
(pIHll with Eq. (IMll . the finite action per unit volume Vd -2 
can be obtained as 



1 = 




A _ 2‘-^q^‘(2s-2d+3) 

(d—1) 2s(d-2)/(2s-l)(2s-d+l) \ 1-2® 

’■+ f 2s-l J 


87r(d — 2) 


(35) 


Having the total finite action, we can use the Brown-York 
method of quasilocal definition m to introduce the fol¬ 
lowing divergence-free stress-energy tensor 

Tab ^ _ T^^ab^ T 2^2(3 J“^ - J7“h 

Stt 

+3a3(5P“'' - P7“'') + 4^7“'' }■ (36) 

Pff 


It is easy to find that d/dt and d/dcj)^ are the Killing 
vectors of the metric m and therefore the associated 
conserved quantities are the mass and angular momentum, 
which are 


M = ^m[{d-l)S^-l], (37) 

Ji = — —Smai. (38) 

iOTT 

We should note that for static solutions {at = 0 and there¬ 
fore S = 1), the angular momentum vanishes, and thus 
the ads are the rotational parameters of the spacetime. 

In order to complete this section, we should calculate 
the entropy of the black branes. Since the area law is ap¬ 
plicable to black holes (branes) in Einstein gravity 

we calculate the entropy through the use 
of Gibbs-Duhem relation 

S = -QU- I. (39) 

It is strightforward to calculate the entropy per unit vol¬ 
ume Vd_ 2 by use of Eq. dMl) with obtained thermodynamic 
and conserved quantities 

5 = (40) 

which confirms the entropy obeys the area law for our case 
(flat horizon). 


Now, we are in a position to check the first law of 
thermodynamics. It has been found that the first law may 
be derived from Smarr(-type) formula for the total energy 
of black holes m- In order to obtain Smarr-type formula, 
one can calculate the mass density as a function of the 
extensive quantities S, J, and Q. Using Eqs. dSOl), (EZl), 
dMl), dlQl) and the fact that r+ is the largest root of /(r), 
we obtain 


M(5,J,Q) = 


[{d-l)Z-l]J 


id-l)l^Z{Z-l) 


(41) 


where J = |J| = and Z = is the positive real 

root of the following equation 


(2s-d+l) 




Z(45) 


d-1 

d-2 


87r(d- 2)J 


Z- 1 


= 0 . 


(42) 


Now, we can regard the mass M{S, J, Q) as a function of 
extensive parameters and obtain the intensive parameters 
conjugate to them in the following manner 


T=(^] 

\9s)j^q 


«7S) 


U=(^-^) 

\dQj 


S,Q \ / s,J 

(43) 

Using the chain rule, it is straightforward to show that the 
intensive quantities calculated by Eq. (031) coincide with 
Eqs. (|2Tp . (j28p sind (EH), and hence we can deduce that 
our black brane solutions satisfy the first law of thermo¬ 
dynamics 


k 

dM = TdS H- '^HidJ, + UdQ. (44) 

5 CLOSING REMARKS 

In this paper, we regarded the first four terms of Lovelock 
gravity with a source of nonlinear electromagnetic field. 
We found a class of rotating solutions which may be in¬ 
terpreted as black branes with inner and outer horizons, 
extreme black branes or naked singularity. 

We used the analytic continuation of the metric to cal¬ 
culate the Hawking temperature and angular velocities of 
the black branes. We then considered the Gauss’s law to 
obtain the finite electric charge of the black branes. In 
addition, we calculated the finite action, mass and angu¬ 
lar momentum by using the counterterm method. Since 
the area law can be applied to black holes (branes) in 
Einstein gravity, we calculated the entropy through the 
use of Gibbs-Duhem relation and found that the entropy 
obeys the area law. 

The Smarr-type formula for the mass was also obtained 
as a function of extensive quantities to compute tempera¬ 
ture, angular velocities and electric potential. It was then 
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confirmed that the thermodynamic and conserved quan¬ 
tities satisfy the first law of thermodynamics. 

Finally, it is also desirable to study the causal struc¬ 
ture, the ratio of shear viscosity to entropy density, ther¬ 
modynamic stability and dynamical properties of the black 
brane solutions derived here. Generalization of these solu¬ 
tions to various (nontrivial) horizon topologies, remain to 
be carried out in the future. In addition. It would be also 
interesting to investigate the effect of curvature-cubed La- 
grangian [sni^ and quartic quasitopological gravity [sg 
with a nonlinear source. We hope to address these issues 
in the future. 

Acknowledgement: We thank the referee for constructive 
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search Institute for Astronomy & Astrophysics of Maragha 
(RIAAM). 
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